Axiverse cosmology and the energy scale of inflation 
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Ultra-light axions (m a < 10 _18 eV), motivated by string theory, can be a powerful probe of the 
energy scale of inflation. In contrast to heavier axions the isocurvature modes in the ultra-light 
axions can coexist with observable gravitational waves. Here it is shown that large scale structure 
constraints severely limit the parameter space for axion mass, density fraction and isocurvature 
amplitude. It is also shown that radically different CMB observables for the ultra-light axion 
isocurvature mode additionally reduce this space. The results of a new, accurate and efficient 
method to calculate this isocurvature power spectrum are presented, and can be used to constrain 
ultra-light axions and inflation. 
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Introduction- Axions pQ are a leading candidate for the 
dark matter (DM) component of the Universe. Proposed 
to solve the strong CV problem, they are also generic in 
string theory [2 , leading to the idea of an axiverse [3]. 
The number of axions in the axiverse is expected to be 
large, due to the topological complexity of string com- 
pact ificat ions. 

In the axiverse the axion mass is given by m 2 a = A 4 // 2 . 
The scales, f a and A a , are both determined separately 
for each axion, and depend on the action, 5, due to non- 
perturbative physics on the corresponding cycle: f a ~ 
Mpi/S^A^ = it 4 e~ s , where M p i is the reduced Planck 
mass, Mp t = 1/SttG. The hard non-perturbative scale is 
\i and its value should be roughly given by the geometric 
mean of the Planck scale and the SUSY scale. Solving 
the strong CV problem with one of the string axions re- 
quires S > 200 [2, 3 , giving rise to stringy values of 
fa ^ 10 12 GeV, near the GUT scale. The exact value of 
*S, however, scales with the area of the corresponding cy- 
cle (a dynamically distributed quantity in the landscape) , 
so that small variations in the area lead to exponential 
variations in the scale of the potential. Thus the mass 
also depends exponentially on the modulus determining 
5, which is evenly distributed, and so the axion mass 
spectrum can be taken as an flat distribution on a loga- 
rithmic scale [3 . 

An important distinction between QCD axions and 
lighter axions is that the temperature dependence of the 
mass drops out and this changes the scalings between 
misalignment angle and relic density. 

If axions are very light, with mass m a < 10 _18 eV 
(ultra-light axions, or ULAs), coherent oscillations of the 
field lead to the suppression of small scale power, and 
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FIG. 1: Adiabatic matter power spectra, with varying ax- 
ion mass m a = 10" 28 , 10" 26 , 10" 25 , 10" 23 eV at fixed den- 
sity fraction Q a /Qd = 0.5 (dashed), and varying Q a /^d = 
0.1,0.5, 1 at fixed m a = 10 -25 eV (solid). Spectra are calcu- 
lated using the methods of Ref. [14] , 



distinguish axions from cold (C)DM [3HB]. Heuristically, 
this scale is the geometric mean of the axion de Broglie 
wavelength and the Hubble scale. Depending on the ax- 
ion mass, this scale could be macroscopic, thus affecting 
observed CMB anisotropics and galaxy clustering power 
spectra. For the classic QCD axion (m a ^ 10 -6 eV), 
this scale is not cosmologically relevant. In the WKB 
approximation (averaging over the fast time scale associ- 
ated with the axion mass, m" 1 ), the axion may be accu- 
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rately treated as a fluid, with sound speed 



cl = 



1 



if k <C 2m a a, 
if k ^> 2m a a. 



(i) 



The scale of structure suppression begins at the scale 
fc m , defined such that those modes with k > k m had 
sound speed (? a = 1 for some time while they were inside 
the horizon [6]. The effect saturates at the smaller scale 
fcj = \/m a H. Therefore, like massive standard model 
neutrinos or more novel warm (W)DM candidates [7], 
axions exhibit suppressed structure on small scales, as 
shown in Fig. [l] 

Axions in inflationary cosmology carry isocurvature 
fluctuations [8], further distinguishing them from ther- 
mally produced CDM. The amplitude of these fluctua- 
tions is set by the energy scale of inflation and is thus 
tied to the amplitude of primordial gravitational waves. 
Axions thus offer an interesting window into the infla- 
tionary epoch, the string landscape, and the multiverse 
[TT] . For a QCD axion, the isocurvature bounds imply 
that tensor modes are unobservably small: surprisingly 
this does not happen for ultra-light axions, and we will 
shortly explain why. 

In past work, these two aspects of axion cosmology - 
the suppression of clustering and the existence of isocur- 
vature perturbations - have been viewed in isolation. 
Probes for axion-seeded isocurvature have been restricted 
to the QCD axion, for which structure suppression on 
small scales is observationally irrelevant [11] [12] . Mean- 
while, observations of CMB anisotropics and galaxy clus- 
tering place limits to axion-induced structure suppression 
[5 , but do not yet include the isocurvature constraint. 

Isocurvature perturbations, gravitational waves and the 
CMB- It is well known that the tensor-to-scalar ratio, 
r, is a probe of the inflationary energy scale and may 
be measured using CMB B-modes [9]. The isocurva- 
ture amplitude of axions is directly related to r, be- 
cause both ULAs and gravitons are massless during in- 
flation. The standard formulae for the tensor, Vh, an d 
scalar, Vki power give the well known result for the 
tensor-to-scalar ratio r = Vh/Vn = 16e, where e is a 
slow roll parameter. Given that the scalar amplitude, 
A s = (l/2e)(Hj /27tM p i) 2 , is well measured, a measure- 
ment of e constitutes a measurement of the Hubble scale 
during inflation, Hi. The isocurvature fraction also de- 
pends on e: measuring it constrains a function of Hj 
and the axion initial misalignment angle [8 . To mea- 
sure Hi using isocurvature one must either constrain or 
make assumptions about the axion initial misalignment 
angle. Tensor modes only measure Hi if they have an 
inflationary origin. The sensitivity to tensor modes and 
isocurvature improves with results from Planck [TO] . 

Isocurvature perturbations are entropy fluctuations of 
the form Sij = (Srii/ni) — (Srij/rij), where the Sni and 
rii are number density fluctuations and average number 
densities, respectively, in each species present. Entropy 
fluctuations arise if there are (nearly) massless spectator 



fields present during inflation [13 , and the axion is one 
example. 

Since the axion is an independent quantum field from 
the inflaton, energetically subdominant during inflation, 
the axion seeds isocurvature perturbations that are un- 
corrected with the dominant adiabatic fluctuations. The 
axion isocurvature fluctuations generated in this manner 
are unavoidable in any standard inflationary scenario as 
long as neither the inflationary fluctuations of the axion 
nor reheating restore the Peccei-Quinn symmetry [TT] . 
This is often the case for the large, stringy, values of 
/„ > 10 12 GeV. 

The de Sitter space quantum fluctuations of the axion 
field (/>, have magnitude 1 : 
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There are constraints (e.g. WMAP9 [12]) on the rela- 
tive amplitude, a, of CDM isocurvature fluctuations de- 
fined by: 



1 



Vn(ko) 



a 



(3) 



where Vs is the isocurvature primordial power spectrum 
evaluated at pivot wavenumber fco- 

The axion power spectrum is given by: 



(Hj/Mpt) 2 




m a a osc 



(4) 
(5) 



The initial misalignment angle is Oi = fa/ f a \ it is fixed 
by the relic density and a osc , which is a function of axion 
mass defined by 3H(a osc ) = m a [6j[T4]. Subsequent to 
a osc the axion redshifts as matter, but displays suppres- 
sion of structure formation. 

Before we can relate a to ra a , Hi and Q a using Eqs. Q- 
([5]), we must clarify the isocurvature normalisation. The 
usual CDM isocurvature normal mode is defined by tak- 
ing S c = 1 as the the initial amplitude of the CDM over- 
density, and normalizing the power spectrum such that 
Vs = V c , where V c is the power spectrum of the CDM 
fractional overdensity. 

If axions are now included as a sub-component of the 
DM with the same equation of state and sound speed as 
CDM (as in Ref. [12] and others), then there is a single 
DM effective fluid (with fractional density perturbation 
Sd) whose isocurvature normal mode is defined by Sd = 
1. If axions carry isocurvature fluctuations, while the 
CDM itself carries only adiabatic fluctuations, then Vs = 



1 These fluctuations also set the variance on initial misalignment 
angle and may alter the axion abundance Q a |H| - In our mass 
range of interest this effect is negligible. 
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FIG. 2: Projected constraints in the {ra a , Q a /^d} plane. The 
region above the black solid lines, labeled "Ruled out 95% 
C.L. LSS (2005)", uses the 2a constraints of Ref. 5 ]. The 
shaded regions are between dashed contours having {0.01 < 
r < 0.1,0.01 < acDM < 0.047}, so that both are potentially 
observable, with the darker shading also allowed by Ref. [5]- 
In these regions it is possible to unambiguously infer Hi from 
a combination of tensor and isocurvature measurements in the 
CMB, combined with a LSS measurement of Q a -The dashed 
black line ("CMB") estimates the modified {a,r} contours 
taking into account isocurvature power suppression for low 
masses (see Fig. |3j). 

(Q a /Qd) 2 T^a, where T a is the axion perturbation power 
spectrum. In the treatment we develop in Ref. [14], we 
incorporate ULAs as a separate effective fluid component 
with their own independent equation of state and sound 
speed, in addition to the CDM. In the axion isocurvature 
normal mode, the initial fractional axion over density is 
8 a = 1, giving Vs = V a - This yields two definitions of 
a, which we call q^cdm (if axions are just included in the 
overall CDM density) and a a (if axions are treated as 
a separate species). The WMAP 9-year constraints to 
axions [12] are derived and stated in terms of q^cdm- 

The two different definitions for the isocurvature frac- 
tion are given by 

a a _ 8e / ttd\ 2 q;cdm , fi x 

i - a a ~ ^i/M vl y ~\n a J i - a C DM ' 1 ) 

Measuring the set {a, A 3 , O a , m a } allows one to constrain 
Hj/Mpi. For any definition of a one has the well defined 
prior range a £ [0, 1]. 

In axion isocurvature models, once m a and Q a are 
specified, r (and thus e and Hi) is uniquely determined 
by a, and vice versa. We visualise the interplay of tensor 
and isocurvature constraints through a schematic plot, 



Fig. [2j which plots contours in r across the entire range 
of cosmologically relevant ra a , using Eqs. ([5|, ([6| for a 
given Q!cdm to fix e at each point. In particular, when 
&osc > a e q 7 where a eq is the scale factor at matter radia- 
tion equality, one finds that r no longer depends on m a at 
fixed a, and so constraints from ULAs can be markedly 
different from CDM axions. The two dashed lines span 
the observable range for ckcdm and r. The isocurvature 
range is 0.01 < c^cdm < 0.047 where the upper bound is 
from Ref. [12 and the lower bound is the forecasted sen- 
sitivity of a cosmic variance limited all sky CMB experi- 
ment in temperature and polarisation [15] . The range for 
tensors is 0.01 < r < 0.1 implying a range of sensitivity 
of an order of magnitude to the energy scale of inflation. 

Fig. [2] shows contours of fixed r and a. Areas shaded 
between these contours have both of observable magni- 
tude. This is in contrast to the QCD axion, due to the 
different scaling of the relic density, and the very low 
mass. In the regions of the {m a ,Q a /Qd} not shaded by 
our contours for r and a, there are two possibilities: ei- 
ther r or a must be unobservable. If high mass ULAs 
exist and constitute a sub-leading fraction of the dark 
matter, bounds to a imply unobservable r, and probe 
low-scale inflation [TT] . Novel to the case of ULAs with 
&osc > tt e q 7 however, is the fact that if they exist and 
are energetically important today, existing bounds to the 
tensor amplitude imply unobservably small a. The oppo- 
site behaviour comes from the switch in the dependence 
of the relic density on mass at a osc = a eq . 

Fig. [2] also shows the constraints to {ra a , Cl a /Cld} from 
Large Scale Structure (LSS) taken from Ref. [5 . Areas of 
the {ra a , Vta/Vtd} plane below the contours of Ref. [5] are 
permitted. These constraints severely limit the region 
where both r and a are simultaneously observable. 

The dark shaded regions in Fig. [2] are particularly in- 
teresting; both regions correspond to simultaneously ob- 
servable values of a and r, while also being consistent 
with the constraints to Q a /Qd of Ref. [5 . With improv- 
ing data that may allow for a detection of Q a [6], it will 
be possible to break the usual degeneracy in {Hi,Q a } 
from constraints to a. In the shaded regions one can 
use an ^ a detection to predict an observable a from an 
observed r and vice versa, thus providing a non-trivial 
cross-check on the inflationary origin of these modes, and 
thus on Hi. Given that there are sources of observable 
tensor modes possible even with low-scale inflation [16] 
these regions provide a novel and truly unambiguous way 
to measure the energy scale of inflation using the con- 
cordance of {a,r, Sl a }: we will present constraints in a 
forthcoming paper [T4] . 

So far we have assumed that constraints to q^cdm will 
map over to constraints to a a . For adiabatic fluctuations, 
the effect of subdominant axions on the CMB observ- 
ables is very small. For isocurvature fluctuations, how- 
ever, the radically different super-horizon solutions [14] 
of axion isocurvature lead to sharply different behavior 
from the more familiar pure CDM isocurvature. This 
mode, as well as the more general suppression of small- 
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FIG. 3: CMB axion isocurvature power spectrum, with adi- 
abatic ACDM for scale (black dashed). We demonstrate the 
normalisation difference between ocdm (grey dot-dash) and 
a a (solid), with Q a /^d = 0.01 implying a normalisation dif- 
ference of (0.01) 2 = 10 -4 . We also show small-scale power 
suppression by the lightest axions. The axion masses are 



scale structure in ULA models, is carefully implemented 
using a modified version of CAMB [17] and is described in 
Ref. [14]. In this case, all other species fall into the grav- 
itational potential wells set up by axions, and so axions 
drive the behavior of the observables, leading to far more 
dramatic effects. We show example spectra in Fig. [3J 

Fig. [3] demonstrates that in the isocurvature mode, 
CMB power is suppressed on small scales (large t\ with 
the scale of power suppression becoming larger as the 
axion mass decreases, just as in P(k) (c.f. Fig. [if. As 
the axion mass increases the axion isocurvature spectra 
asymptote to CDM-like behaviour. 

The suppression of power will be important for ULAs 
in altering the isocurvature constraints. Since the isocur- 
vature power spectrum falls off rapidly at large i, most 
constraining power on isocurvature comes from the ad- 
dition of power along the low-£ plateau before the first 
peak at i ~ 200. When the isocurvature power is sup- 
pressed along this plateau the isocurvature spectrum re- 



mains significant only at lower and lower £. Therefore we 
should expect that not only will allowed values of a a be 
different from c^cdm due to normalisation, but also due 
to the power suppressing properties of ULAs. The effect 
of this is estimated from the reduced number of modes 
available to measure isocurvature fraction and is shown 
in Fig. [2] Isocurvature becomes harder to measure and 
further constrains the observable region for {a, r} at the 
lowest masses, m a < 10 -28 eV. The lowest mass region is 
harder to access observationally using LSS measurements 
since the structure suppressing properties of the axions 
only occur on very large scales [6j. In addition, producing 
an observable relic density with m a < 10 -28 eV would re- 
quire additional physics: for example a large number of 
axions with nearly degenerate masses. 

Conclusions- In this letter we have demonstrated that 
in the case of ultra-light axions one is able to unambigu- 
ously infer the energy scale of inflation from their isocur- 
vature fraction by using large scale structure constraints 
to bound the relic density. In addition, there are regions 
of parameter space allowed by current constraints where 
both the isocurvature fraction and the tensor to scalar 
ratio are within observable reach of near future CMB ex- 
periments. This predicted concordance of three observ- 
ables is a potentially powerful probe of the energy scale 
of inflation. In the context of the axiverse, the inferred 
value of Hj from observed tensor modes would predict 
observable axion isocurvature across more than four or- 
ders of magnitude in axion mass. We present constraints 
to this model in a forthcoming paper [14]. 
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